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Abstract — Constellation Constrained (CC) capacity regions 
of two-user Gaussian Multiple Access Channels (GMAC) have 
been recently reported, wherein an appropriate angle of rotation 
between the constellations of the two users is shown to enlarge 
the CC capacity region. We refer to such a scheme as the 
Constellation Rotation (CR) scheme. In this paper, we propose 
a novel scheme called the Constellation Power Allocation (CPA) 
scheme, wherein the instantaneous transmit power of the two 
users are varied by maintaining their average power constraints. 
We show that the CPA scheme offers CC sum capacities equal (at 
low SNR values) or close (at high SNR values) to those offered 
by the CR scheme with reduced decoding complexity for QAM 
constellations. We study the robustness of the CPA scheme for 
random phase offsets in the channel and unequal average power 
constraints for the two users. With random phase offsets in the 
channel, we show that the CC sum capacity offered by the CPA 
scheme is more than the CR scheme at high SNR values. With 
unequal average power constraints, we show that the CPA scheme 
provides maximum gain when the power levels are close, and the 
advantage diminishes with the increase in the power difference. 

Index Terms — Constellation constrained capacity regions, mul- 
tiple access channels, power allocation, finite constellations. 



I. Introduction AND Preliminaries 

In traditional networks, information exchange between the 
subscriber devices and the base station is realized by schedul- 
ing the transmissions in disjoint set of point-point channels, 
i.e., by separating the subscribers through TDMA, FDMA, 
or CDMA. Over the last few decades, enormous research has 
taken place to look beyond cliannel separation, and develop 
advanced physical layer techniques that can provide capacity 
gains at the cost of joint-processing of the signals from 
different subscribers [l]-[7]. Along that direction, determin- 
ing the capacity regions of multi-terminal networks [8], and 
designing appropriate low-complexity signalling schemes [9] 
has received a lot of attention. 

Till date, capacity regions are known only for a certain class 
of network configurations such as the multiple access channels 
(MAC) [3], [4], and Gaussian broadcast channels (BC) [10], 
to name a few. Though, contributions on such configurations 
started for channels with few users (such as two-user MAC/BC 
or three-user MAC/BC), the results have now been generalized 
to arbitrary number of users [6], [7], [11]. Further, for Gaussian 
MAC (GMAC), the capacity achieving input distribution is 
known to be Gaussian, which is a continuous distribution [8] . 
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However, in practice, the input constellations are finite in size 
and are uniformly distributed. As a result, the known results for 
continuous input does not shed light on the actual performance 
of finite constellations. This difference in the nature of the 
input alphabet has motivated researchers to revisit the GMAC 
(and similar channels), and study them from the view point of 
finite input constellations. As a first step, some results have 
been reported for channels with few users [13]- [27] . This paper 
is along that direction, and we deal with two-user GMAC with 
finite input constellations. 

Influenced by some preliminary works on MAC with finite 
inputs constellations [13], [14], a detailed study on the CC 
capacity [28] regions of two-user GMAC has been reported in 
[15]. It is shown in [15] that introducing an appropriate angle 
of rotation between the constellations provides enlargement 
in CC capacity region. Further, angles of rotation which 
maximizes the CC sum capacity have been provided for 
some known constellations. We refer to such a method of 
enlarging the CC capacity region as the Constellation Rotation 
(CR) scheme. Note that the CR scheme is a Non-Orthogonal 
Multiple Access (NO-MA) scheme, wherein the two users 
transmit during the same time and in the same bandwidth. 
It is also shown in [15] that the CC capacity region of the 
CR scheme strictly encloses the CC capacity region of the 
FDMA and the TDMA. It is highUghted in [15] that the 
above behaviour is not observed with Gaussian inputs, which 
in turn shows the importance of studying these channels with 
finite input constellations. Other than two-user GMAC, similar 
research on finite input constellations have also been reported 
for MAC with quantization and fading [16]-[18], Gaussian 
broadcast channels [19], interference channels [20], [21], relay 
channels [22], [24], point-point MIMO channels [26], and 
secrecy channels [27]. 

In the schemes proposed in [15]-[27], either the channels 
are assumed to be fixed, or some form of channel state 
information is available at the transmitters. For instance, the 
CR scheme [15] assumes fixed channels, and the technique 
is sensitive to the choice of the relative angle of rotation. 
If the channels introduce random phase offsets (say, due to 
clock synchronization problems), then the resultant relative 
angle need not provide maximum gains. In such a case, even 
if the random phase offsets are made available to the receiver, 
the CR scheme becomes ineffective as the transmitters do not 
have the knowledge of the phase values. As a result, there 
is a need for signal design which is robust to random phase 
offsets. In this paper, we propose a NO-MA scheme called the 
Constellation Power Allocation (CPA) scheme [16], wherein 
instead of introducing rotation between the constellations, we 
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vary the transmit powers of the two users. Unlike the CR 
scheme, the proposed CPA scheme is robust to random phase 
offsets in the channel. The contributions of this paper are 
summarized as below: 

• We propose a novel transmitter side technique called the 
CPA scheme to obtain enlargement in the CC capacity 
region of two-user GMAC. In the proposed scheme, the 
transmit power of the two users are varied while retaining 
the average power constraint for each user 

• For the CPA scheme, the transmit power of the users are 
varied through a scale factor a e [0, 1]. For such a model, 
we propose the problem of finding an appropriate a that 
maximizes the CC sum capacity. For simpler computa- 
tions, we propose a deterministic metric to compute the 
scale factor such that the CC sum capacity is maximized 
at high Signal-to-Noise Ratio (SNR) values. We compute 
the CC sum capacities for some known constellations, 
and show that the CC sum capacities offered by the CPA 
scheme are equal (at low SNR values) or close (at high 
SNR values) to those offered by the CR scheme in [15]. 
(Section Imb 

• For regular QAM constellations, we first identify that the 
CPA scheme provides a sum constellation Ssum whose 
in-phase and the quadrature components are separable, 
and subsequently show that the CPA scheme provides 
lower decoding complexity than the CR scheme. This 
advantage is shown to come with no significant reduction 
in the CC sum capacity. To exploit the reduced decoding 
complexity offered by the CPA scheme, we propose 
independent coding along the in-phase and the quadrature 
components for each user We also propose TCM based 
code pairs to approach the CC sum capacity for 16- 
QAM constellation. We point out that the low decoding 
complexity advantages of the CPA are applicable only for 
GMAC with no random phase offsets. (Section llVb 

• We study the robustness of the CPA scheme for two-user 
GMAC with random phase offsets in the channel. For 
such channels, it is clear that the CR scheme does not 
improve the CC sum capacity due to the random phase 
offsets unknown to the transmitters. We show that the 
CC sum capacity offered by the CPA scheme is more 
than the CR scheme at high SNR values. We also study 
the robustness of the CPA scheme for unequal average 
power constraints for the two users, where we show that 
the gains of the CPA scheme decreases as the power 
difference between the users increase. (Section |V) 

Notations: Throughout the paper, boldface letters and capi- 
tal boldface letters are used to represent vectors and matrices, 
respectively. For a random variable X which takes value from 
the set S, we assume some ordering of its elements and use 
X{i) to represent the z-th element of S, i.e., X{i) represents 
a realization of the random variable X. We use the symbol z 
to represent \/— 1. Cardinality of a set S is denoted by \S\. 
Absolute value of a complex number x is denoted by 
and E [x] denotes the expectation of the random variable x. 
A circularly symmetric complex Gaussian random vector, x 
with mean fi and covariance matrix F is denoted by x 



cscg (/x,r). 

II. Two-user GMAC: Signal Model and CC Sum 
Capacity 

The model of two-user GMAC consists of two users that 
need to convey independent information to a single destination. 
It is assumed that User-1 and User-2 communicate to the 
destination at the same time and in the same frequency band 
(the two users employ a NO-MA scheme). Symbol level 
synchronization is assumed at the destination. In this section, 
we assume no random phase offsets introduced by the channel. 
The two users are equipped with finite complex constellations 
Si and ^2 of size iVi and N2, respectively such that for 
Xi G Si, we have = 1. Let Pi denote the average 

power constraint for User-i. When User-1 and User-2 transmit 
symbols \fP\x\ and \fP2x1 simultaneously, the destination 
receives the complex symbol y given by 

y = y/FiXi + y^X2 + z, where z - CSCG (O, a^) , (1) 
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and ^ is the variance of the AWGN in each dimen- 
sion. The CC sum capacity of two-user GMAC [15] is 
I {VPiXi + \fP2X2 '■ y) , which is given in (O at the top of 
the next page. 

In [15], an appropriate angle of rotation between the con- 
stellations is shown to increase the CC sum capacity. In this 
paper, we introduce the CPA scheme to increase the CC sum 
capacity. Before introducing the CPA scheme, we recall the 
definition of the sum constellation and uniquely decodable 
(UD) constellation pairs [15]. Given two constellations Si and 
^2, the sum constellation of Si and 52 is given by 

Ssum — \^\fPlXl + \/~P2X2 I V Xi € Si,X2 £ ^2 1 . 

The adder channel in the two-user GMAC can be viewed as 
a mapping given by 

(j) {yKxi, ^/P2X2^ = ^/Kxi + ^/P2X2. 

A constellation pair (iSi, ^2) is said to be uniquely decodable 
(UD) if the mapping 4> is one-one. 

III. Constellation Power Allocation Scheme for 
Two-User GMAC 

In Section Ull we have assumed unit average power on Si 
and 52. Therefore, to meet the average power constraint Pi, 
User-i can transmit symbols of the form \/P'iXi. We now 
propose an alternate method to transmit the symbols of Si 
by maintaining the average power constraint Pi. To explain 
the new scheme, we let T = {1, 2, 3, • • • , T — 1, T} denote 
the set of the indices of complex channel use, where T is 
the total number of channel uses. Assuming T to be an 
even number, let T^^j^j = {1, 3, 5 • • • , T — 3, T — 1} and 
Tiven = {2,4,6- • • ,r - 2,r} denote the set of odd and 
even indices of channel use, respectively. We use the variable 
t to denote the instantaneous channel use index. We also use 
a € [0, 1] to denote a real valued variable which is used to 
vary the transmit power of each user Using the pseudo code 
representation, we explain the CPA scheme below. 
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JVi-l iV2-l 
fel=0 fc2=0 



log2 



I{^/pIxl + ^/P^X2 : y) = \og2(NiN2)- 
exp(-|z|2/^2) 



(2) 



iVi ~ 1 iV2 - 1 

= E E i°g2 

fcl=0 fc2=0 



iVi-1 ]V2-1 

^ 1^ exp(-|y(2^^S)Pr(a:i(fei)-xi(n)) + Vs^(x2{fc2)-a;2(i2))|V2<T^ 

il=0 12=0 



(3) 



-f'^Hv^a;! + ^1^X2 : y) ■ 



iVi-l JV2-I 

E E ^ 



log 



iVi-1 JV2-I 

1^ 1^ exp(-|yPZa-i{fci) + yp;x-2(fc2)- yPZ'^^iCn)- V^a;2(i2) + 2|Vf^^ 
n=o 12=0 

V 

/3(fcl,fe2,z) 



(4) 



A(fei,fe2) 



IF t e Todd 

• User-l transmits ^(2 — a)Pia;i for xi £ 5i 

• User-2 transmits ^/aP2X2 for 0:2 G ^2 
ELSE 

• User-l transmits yJaPiXi for a;i G 5i 

• User-2 transmits ■\/(2 — a)P2a;2 for X2 G ^2 
END 

During odd indices, the destination receives a symbol yt 
of the form yt — ■\/(2 — q;)Pi.ti + -y/(x^a:;2 + z. Similarly, 
during even indices, the destination receives a symbol yt of 
the form yt = ^/aPixi + ^(2 — a)P2X2 + z. With the above 
mentioned power allocation, the average power for User-l is 
(2 — a)Pi and aPi in odd and even channel uses, respectively. 
Similarly, the average power for User-2 is aP2 and (2 — a)P2 
in odd and even channel uses, respectively. With this, the 
average power Pi is maintained for User-i. When the two users 
employ identical constellations and equal average power, an 
appropriate value of a can provide the UD property at the 
receiver for every channel use. 

We now proceed to find the optimal a that maximizes the 
CC sum capacity of the CPA scheme. 

A. Optimal a for the CPA Scheme 

We consider the identical constellation case, i.e.. Si = 
S2, and focus on finding a such that the CC sum 
capacity is maximized. Through the CPA scheme, the 
two users switch the scale factors on alternate chan- 
nel uses. During odd channel uses, the destination views 
the sum constellation 5s„,„,odd given by Ssu7n,odd - 
|-\/(2 — a)PiXi + \foLP2X2 I Va;i,a;2| • During even chan- 
nel uses, the destination views Ssum.eve.n given by 
Ssura.evev, = ^VaP[xi + ^(2 - a)P2a;2 I Va;i,X2| • If the 
two users have equal average power constraint, the destination 
views the same sum constellation on every channel use. On 
the other hand, if the two users have unequal average power 
constraint, then the sum constellation seen by the destination 
is not the same during the even and odd indices of the channel 



use. Considering the general case of unequal average power 
constraints, the CC sum capacity of the CPA scheme is given 
by ^ 

- ^ / (v/(2 - (^)Pixi + V^2X2 ■■ y) , 

where il = {a, 2 — a}. Note that the CC sum capacity can 
be increased by choosing an appropriate a e [0, 1]. However, 
we know that a = corresponds to single-user transmission. 
As a result, henceforth, we consider selecting an appropriate 
a in the interval (0, 1]. Thus, the CC sum capacity can be 
maximized by choosing aopt given by, 

"opt = argmax 1:^1 (\/ (2 - a)PiXi + ^ aP2X2 : . 

(5) 

Note that the above objective function involves expectation of 
a non-linear function of the random variable z, and hence, its 
closed form expression is not available. Therefore, computing 
Oopt is not straightforward. For high values of ^ and the 
following theorem provides a deterministic metric (which is 
independent of the variable z) to choose a such that the CC 
capacity is maximized. Henceforth, high SNR values imply 
high values of ^ and 

Theorem 1: At high SNR values, the optimum scale factor 
a required to maximize the CC sum capacity is approximated 
closely by a* where 

a* = argmin Q(q:), 
"e(o,il sen 

where Q{oi) is given in ^ at the top of this page. 

Proof: Since A^i and N2 are constants, we have the 
following equality, 

argmax ^ (2 — a)Pixi + \fdLP2X2 ■ y) 
"e(04] ago 

= argmin (2 - a)Pixi + \J 6LP2X2 : y), 

"£{0,1] ago 

where /'^^ (-^Z (2 — (x)P\Xx + \/^P^X2 : y) is given in (|4|, and 
we use Pl to denote (2 — a)Pi, and P5 to denote aP2. Note 
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that the individual terms A(A:i , fc2) of /'^^ {\fPLXi + \fPsX2 ■ 
y) are of the form E [log2(/3(fci, ^2, z))] for a random variable 
/3(fci, fc2, z). Applying Jensen's inequality: 

i?[log2(/3(fci, fc2, Z))] < l0g2[i?(/3(fci, k2,z))] 

on A(fci,fc2), we have 

/(I) ( v/P^xi + /f^a;2 : y) < Q(a), (6) 

where (5(q;) is given by (|3]l. In the rest of the proof, we show 
that at high SNR values, the following approximation holds 
good: 

argmin Q{a) « argmin /^^^(v^P^xi + \/PsX2 : y). 

Note that the term I'-^^y/PLXi + \fPsX2 '■ y) can be written 
as in ^ and ((SJ where 

/X(fci,fc2,«l,«2) =V^2:i(fci) + \/P^X2{k2) 

- \/PLXi{ii) - \/PsX2ii2), 

and M(fci,fc2) = |Al(fci,^2)| such that Al(A;i,fc2) is given 
by 

M{ki, k2) = {(ii, 12) 7^ (^1,^2) I /^(ii, 12, fci, ^2) = 0}. 

Removing independent terms of the form exp ^— ■'^^ in ([S]), 
we have 

argmin l'^^^{\fi^xi + \fPsX2 : y) 

= argmin /'^^^(yP^xi + \fPsX2 ■ y), 

"e(o,i] 

where J*-^-' [y/ PlXi + y/PsX2 : y) is given in ( fTOl i. Further, at 
high SNR values, we have the approximation, 

argmin V I'''^\^/pIxi + y/P^X2 : y) 

w argmin ^ /^^^ (a/pJ^X'i + \/PsX2 : y), 

"e(o,i] sen 

where I*^^^ (y/PLXi + y/PsX2 '■ y) is given by (fTTI ). At high 
SNR values, each term 7(^1,^2,2) in (fTTl i is small, and 
hence we use the approximation log2(l + 7(fci,fc2,z)) « 
log2(e)(7(fci, k2, z)) to obtain (fTZt . Evaluating the expectation 
in ( fT2] i. we get ( fTsT l. Once again, applying the approximation 
log2(l + S{ki, k2)) ~ log2(e)(5(fci, ^2)) in we get (O, 
which is denoted by J'^^-' (VPf^a^i + \rPsX2 : y). 

Now, we consider the term (V"PLa;i + v^^.T2 : y) given 
in (fTSI l and show the following equality: 

argmin I^^\\/PlXi + \/PsX2 : y) 

= argmin V" /'^'^^ (^P^a;! + V^2;2 : y). (9) 

Once the above equality is proved, the statement of this 
theorem also gets proved since I^^^ (V P^xi + y/PsX2 : y) is 
a scaled version of Q{a) (as shown in (fTSIl). Towards proving 
the equality in (|9]l, note that at high SNR values, 6{ki, ^2) is 
small for all values of a. For those values of a which provide 



the UD property, we have log2(l + M{ki, k2)) = Vfci, k2. 
However, for those values of a which do not provide the 
UD property, we have log2(l + M{ki,k2)) > 1 for some 
fci, fc2- Further, at high SNR, for all a, log2(l + 5(fci, fc2)) << 
1, Vfci,fc2- Due to these reasons, the values of a which do 
not provide the UD property do not minimize I^^\-\/7lXi + 
\fPsX2 ■ y) as well as I^'^\\/Pi,xi + \fPsX2 ■ y)- As a 
result, the optimal value of a must belong to the set which 
provides the UD property. For such values of a, we have 
I^^HVPExi + ^X2 : y) = /^(/P^xi + V^X2 : y), 
and hence the equality in (|9]l holds. This completes the proof. 

■ 

Using the results of Theorem [T] for high SNR values, we 
propose to find a* which minimizes X^sgO Qi'^)' ^ '^ight 
upper bound on '^g^en ^'^'(V^^t^i + \fPsX2 '■ y)- However, 
note that for small to moderate values of SNR, the values of a* 
obtained by solving ^ need not maximize X^gg^ ^iV^LXi + 
VPsX2 '■ y) since the bound in ^ is not known to be tight. 
It is also clear that solving Q is easier than solving (|5) since 
Saeo Qi*^) is deterministic and independent of the term z. 

B. Numerical Results 

In this section, we first compute a* for the case of equal 
average power constraint for the two users (i.e.. Pi = Pj)- 
For the simulation results, we use cr^ ~ 1 and SNR = Pi. 
The values of a* are obtained by varying a from to 1 in 
steps of 0.01. In Table [H the values of a* are presented for 
some known constellations. The CC sum capacities of QPSK 
and 8-PSK are also provided in Fig. [T| and Fig.|2] respectively 
for the following schemes: (/) CPA scheme, (//) CR scheme, 
{Hi) neither CPA nor CR, and (/v) with both CPA and CR. For 
the scheme "with both CPA and CR", the pairs {a* , 6*) are 
computed using a metric which is obtained on the similar lines 
of Theorem [T] From the figures, note that the CPA scheme 
provides CC sum capacities close (or equal) to the CR scheme 
for all SNR values. Note that the scheme "with both CPA and 
CR" does not provide significant CC capacity gain over the 
CPA scheme. 

QPSK Constellation 




5 10 15 20 25 30 

SNR in dB 



Fig. 1. CC sum capacity of QPSK constellation with equal average power 
constraint and no random pliase offsets. 

To verify the results of Theorem [T] we plot the CC sum 
capacity of QPSK and 8-PSK with aopt and a* in Fig. 
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fel=0 fe2=0 



log2 



JVi-1 N2-I 
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fcl=0 fc2=0 



log2 



Afi-l iVa-l 

exp(-|z|2/a2)(l + A/(fci,fe2))+ X] E exp(-|M(fei,fc2,n,i2) + 2|V'^^) 
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(V) 
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iVi-l JV2-I 

E E^ 

fcl=0 fc2=0 



log2 (1 +M(fci,fc2)) +log2 



iVi-1 JV2-I 



1 



(1 + M(fcl,fc2)) 



E E 



il=0 12=0 
/j(fcl,fe2, 21,12)7^0 



-|M(fci,fc2,ii,i2) + + 1^1^ 



(10) 



ATl-l iV2-l 

E E ^ 

fei =0 fe2=0 



log2(l + Af(fci,fc2))+log2 



1 



(l + A/(fci,fe2)) 



iVi-1 ATj-l 

E E 

il=0 12=0 
H{fcl,fc2,n.»2)7^0 



exp 



fe2, ii, 82) + zf 



(11) 



7(fcl,fc2,z) 




Fig. 2. CC sum capacity of 8-PSK constellation with equal average power 
constraint and no random phase offsets. 



Fig. 3. CC sum capacity of QPSK and 8-PSK constellations with agpt 
a* for GMAC with equal average power constraint and no random phase 
offsets. 



[3] From the figure, we see that a* provides approximately 
same CC sum capacity values not only at high SNR but 
also at moderate SNR values. In addition, we also present 
the corresponding values of Oopt and a* in Table. [Ill which 
shows that aopt and are different for low SNR values, and 
are approximately same for high SNR values. These results 
demonstrate that the metric X^aen Qi'^) is a tight upper bound 
on Eaesi + VPsX2 : y) at high SNR values. 

For the unequal average power case, we compute a* for the 
QPSK constellation. We use SNR = Pi (with = 1). For the 



simulation results, we have considered the following relations 
between Pi and P2: (i) P2 = 0.3Pi, (//) P2 = 0.5Pi and (/;;) P2 
= O.75P1, and (iv) P2 = 0.9Pi. In Fig.g] the CC sum capacities 
of QPSK are provided for the CPA scheme, the CR scheme, 
and the "neither CPA nor CR" scheme. When compared to 
the "neither CPA nor CR" scheme, the CPA scheme provides 
increased CC sum capacities at moderate to high SNR values, 
when Pi and P2 are close. When compared to the CR scheme, 
the CPA scheme provides marginal improvement in the CC 



6 



I^'^H^/P^Xl + ^/P^X2 : y) 



Ni-1 N2-I 

E E 



log2(l + M(fci,fc2))+£; 



log2(e) 



(l + M(fci,A;2)) 



JVi ~i N2-1 

E E 

!1=0 22=0 



exp 



fc2,il,«2) + 



iVi - 1 JV2 - 1 

^ E E 

fcl=0 fc2=0 



log2(l + M(fci,fc2)) + 



1 



(l + A/(fcl,fc2)) 



2^ exp 



il=0 12=0 
M(fcl.'=2.»l.»2)#0 



2^2 



5(fcl.fc2) 



JVi - 1 JV2 - 1 

E E 

fel=0 fc2=0 



log2(l + M(fcl,fc2)) + -log2 



JVi - 1 JV2 - 1 



exp 



-|^(fci, fc2, »1, t2)|^ 
2^2 



JVi-l JV, -1 

E E 

fel=0 fe2=0 



log2(l + M(fcl,fc2))+log2 



1 



(1 + M(fci,fc2)) 



iVi-1 JV2-I 

E E 

n=o 12=0 

(i(fcl,fc2,*l,!2)7^0 



~|^(fci,fc2,n,»2)|^ 

2^2 



^ JV1-1JV2-1 

2 E E 

fel=0 fc2=0 



iVi-1 JV2-I 

^ J2 (-l/PT^^'iCfci) - /Pr^-'iCn) + {^/P^X2{k2) - ^/P^X2{i2))\^ /2a' 
n=0 22=0 
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Q{a) (16) 



sum capacities only at high SNR values, when Pi and P2 
are close. The figure shows that the gains of the CPA scheme 
diminishes when the power difference is more than 3 dB (as 
shown for the case P2 = 0.5Pi and P2 = 0.3Pi). 

C. Reduced ML Decoding Complexity for QAM Constellations 
with CPA 

In this subsection, we highlight the advantage of using 
CPA over CR for the class of regular QAM constellations. 
For uncoded transmission, when CR is employed for QAM 
constellations, the in-phase and the quadrature components of 
the symbols of Ssum are entangled. However, with CPA, since 
the scale factor a is real valued, the in-phase and quadrature 
components of the symbols in Ssum do not get entangled. 
As a result, Ssum is separable, and can be written as the 
cross product of in-phase and quadrature components of its 
symbols. In particular, if the two users employ square regular 
7\/-QAM constellation, then there are AI points in Ssum along 
the in-phase and the quadrature component, respectively. Since 
Ssum is separable, the destination can decode the in-phase and 
the quadrature components independently. Therefore, the worst 
case ML decoding complexity is 0{AI). However, with CR 
scheme, the worst case ML decoding complexity is 0{M^). 
Therefore, for the class of QAM constellations, CPA provides 



lower decoding complexity with negligible loss in the CC sum 
capacity when compared with CR. 

IV. Channel Coding for QAM Constellations with 

CPA 

In this section, we design code pairs based on TCM (Trellis 
Coded Modulation) [12] to achieve sum-rates close to the CC 
sum-capacity of QAM constellations using the CPA scheme. 
Throughout this section, we assume equal average power 
constraints for the two users with Pi = P2 = P- If Si and ^2 
represent regular 7\/-QAM constellations, then the in-phase 
and the quadrature components of Ssum are respectively of 
the form \/{2 — a)Pxii + \/ aPx2i and \/{2 — a)PxiQ + 
\faPx2Q, where xu £ Sii,XiQ e S,q and Su and S^q 
denote the corresponding v'M-PAM constellations for the z-th 
user along the in-phase and the quadrature dimension, respec- 
tively. The set of in-phase and quadrature symbols of Ssum are 
respectively denoted by Ssumj and Ssum,Q, and are given by 
Ssumj = |(V(2 - a)Pxii + VaPx2i) \ xu € , and 

Ssum.Q = {(V(2 - a)PxiQ + VaPx2Q) I x,Q e 5,q| . 
Since the CPA scheme makes the in-phase and quadrature 
components of Ssum separable, the symbols of Ssumj can 
be decoded independent of the symbols of Ssum.Q, thereby 
reducing the decoding complexity. To facilitate this, User-i 
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TABLE I 

Numerically computed a* with equal average power constraint and no random phase offsets. 



SNR in dB 


QPSK 


8-PSK 


16-PSK 


16-QAM 





0.74 


0.65 


0.65 


0.48 


2 


0.65 


0.85 


0.85 


1.00 


4 


0.52 


0.74 


0.74 


1.00 


6 


0.46 


0.67 


0.67 


1.00 


8 


0.43 


0.65 


0.66 


0.84 


10 


0.41 


0.64 


0.67 


0.74 


12 


0.41 


0.59 


0.70 


0.68 


14 


0.40 


0.53 


0.74 


0.46 


16 


0.40 


0.50 


0.79 


0.62 


18 


0.40 


0.49 


0.78 


0.13 


20 


0.37 


0.49 


0.59 


0.12 


22 


0.24 


0.49 


0.58 


0.12 


24 


0.15 


0.49 


0.56 


0.12 


26 


0.10 


0.49 


0.55 


0.12 


28 


0.06 


0.48 


0.55 


0.12 


30 


0.04 


0.13 


0.55 


0.12 



TABLE 11 

Numerically computed values of a with equal average power constraint and no random phase offsets . 



SNR in dB 


"opt QPSK 


a* for QPSK 


"opt 8-PSK 


a* for 8-PSK 





0.88 


0.74 


0.92 


0.65 


4 


0.55 


0.52 


0.83 


0.74 


8 


0.44 


0.43 


0.66 


0.65 


12 


0.41 


0.41 


0.60 


0.59 


16 


0.41 


0.40 


0.52 


0.50 


20 


0.38 


0.37 


0.49 


0.49 


24 


0.19 


0.15 


0.49 


0.49 




10 20 3D 10 20 30 

SNRindB SNR in dB 



Fig. 4. CC sum capacities of QPSK with unequal average power constraint 
and no random phase offsets. 



should make sure that the symbols of SiQ and Su are not 
coded jointly. Therefore, each user can have two encoders one 
along each dimension. Let the subscript X £ {J, Q} denote 
either the in-phase dimension or quadrature dimension. For 
each i £ {1,2}, let User-i be equipped with a convolutional 
encoder Cix with rriix input bits and rriix + 1 output 
bits. Throughout the section, we consider convolutional codes 
which add only 1-bit redundancy. Let the rriix + 1 output 
bits of Cix take values from \/M-PAM constellation Six 



such that \Six\ = 2'"'"^+^. Henceforth, the set of codewords 
generated from Cix and C2X are represented by trellises Tix 
and T2X respectively. 

We assume that the destination performs joint decoding of 
the in-phase symbols of User-1 and User-2 by decoding for 
a sequence over Ssumj on the sum trellis, Tsum.i (see [15] 
for the definition of the sum trellis). Similarly, joint decoding 
of the quadrature symbols of User-1 and User-2 is performed 
by decoding for a sequence over Ssum.Q on the sum trellis, 
Tsum,Q- Due to the existence of an equivalent AWGN channel 
in the GMAC set-up, the sum treUis, Tsum,x has to be labeled 
with the elements of Ssum.x satisfying the design rules in 
[12]. However, such a labeling rule can be obtained on Tsum,x 
only through the pairs (Tix, T2x) and {Six,S2x)- Hence, we 
propose labeling rules on Tix and T2X using Six and S2X 
respectively such that Tsum,x is labeled with the elements of 
Ssum,x as per Ungerboeck rules. 

Since the number of input bits to Cix is rriix, there are 
2mix edges diverging from (or converging to; henceforth, we 
only refer to diverging edges) each state of Tix- Also, as 
there is only one bit redundancy added by the encoder, and 
as \Six\ = 2'"'-^+^, the edges diverging from each state have 
to be labeled with the elements of a subset of Six of size 
2'nix Therefore, for each i. Six has to be partitioned into 
two sets S]-^ and S^^ , and the diverging edges from each state 
of Tix have to be labeled with the elements of either Sj^ or 
S^-^. From the definition of sum trellis, there are 2"'^^+™^'^ 
edges diverging from each state of Tsum,x and these edges 
get labeled with the elements of one of the following sets, 
A = lSl^+S2j^ I i,je{l,2}>. To satisfy Ungerboeck 
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design rules, the transitions originating from the same state 
of Tsum.x must be assigned symbols that are separated by 
largest minimum distance. Therefore, the problem addressed 
is to find a partitioning of Six (a/M-PAM constellations) 
into two sets 5^^^ and S]^-^ of equal cardinality such that the 
minimum Euclidean distance (denoted by dmin) of each one 
of the sets in A is maximized. However, since dmin values of 
the sets in A can potentially be different, we find a partitioning 
such that the minimum of the dmin values of the sets in A is 
maximized. 

A. Designing TCM Schemes with 16-QAM Constellation 

The set partitioning problem described above depends on 
the structure of the sum constellation. As a result, the solution 
to the set partitioning problem depends on the choice of a. 
For arbitrary values of M and a, we are unable to solve 
the set partitioning problem due to lack of structure on 
Ssum of two QAM constellations. However, through computer 
simulations, we have found a solution for the above problem 
for 16-QAM constellation. For 16-QAM constellation, Sn and 
SiQ are of the form Sii = ^{-3,-1,1,3} and SiQ = 
-^^{— 3i, —I, I, 3i}. For this set-up, we obtain a two way par- 
tition of Sii and SiQ such that the minimum of the dmin values 
of the sets in A is maximized. In particular, the optimal parti- 
tion is obtained for different P values. For each value of P, the 
optimal partition is obtained by using the corresponding a* as 
given in the last column of Table H] The optimal partitions for 
all values of P are found to be Slj = -^^{—3,l},Sfj = 
^{-^,3^},SlQ = ^{-3,l}and5,^"L 
which is the Ungerboeck partitioning. With this set partition- 
ing, trellis code pairs based on TCM can be designed in order 
to transmit 2 bits for each user (for each user mix — 1 
bit is transmitted along each dimension) using the 16-QAM 
constellation. 

In the preceding sections, we have studied the advantage of 
the CPA scheme for two-user GMAC with no phase offsets in 
the channel. In the next section, we study the robustness of 
the CPA scheme for random phase offsets in the channel. 



along the similar lines of the model in Section|ll] In addition to 
the steps needed to compute the mutual information values in 
Section ini in this case, we have to take expectation of the mu- 
tual information values over di and 62- Hence, the CC sum ca- 
pacity is given by Eg^^g,^ [l {y/PiXi + \fPiX2 : y \ 6*1, 6*2)], 
where / {\fP\X\ + \fP2X2 '■ V \ ^1,^*2) is given in ( fTSl l. 

Note that the CC sum capacity is a function of the distance 
distribution (DD) of the sum constellation Ssum given by 



V xi e Si,x2 e S2^ 



Despite having random phase offsets in the channel, the 
DD of Ssum can be changed by scaling the input symbols 
of one user relative to the other With such a change in 
the DD of Ssum, the CC sum capacity can be increased 
by choosing an appropriate a. Towards that direction, we 
apply the CPA scheme (given in Section HIH) for this channel. 
With the CPA scheme, the CC sum capacity is given by 
Ee,,e2 [Eaenl{V(^^W\xi + V^X2:y \ 01,^2)], 
where = {a, 2 — a}. Since the CC sum capacity is a 
function of a, we have to compute a such that 



Eg 




y 



is maximized, where Pl = {2 — a)Pi and Ps = aP2- In other 
words, we have to solve the following optimization problem. 



ae{0,l] 



Plxi + VPsx2 ■■ y \ 01, 



(21) 



Note that the closed form expression of Eg^,g^\I{y/PLXi + 
\fPsX2 '■ y I 81,92)] is not available. Therefore, computing 
Oopt is not straightforward. On the similar lines of Theorem 
[T] we use a* = arg minQg(o,i] Z^tseo Qp(q^)' ^° obtain the 
appropriate values of a, where Qp{a) is given in ( fT9] l, and 
/ii(fci, k2,ii, 12) is given in ( |20] |. Note that Qp{a) is free from 
the random variable z. Hence, it is computationally easier to 
obtain a* . 



v. cpa scheme for two-user gmac with random 
Phase-Offsets 

In this section, we consider a two-user GMAC with random 
phase offsets introduced by the channel for both the users 
[14]. Similar to the signal model in Section the two users 
are equipped with constellations Si and ^2 of size A^i and 
N2, respectively. When User-1 and User-2 transmit symbols 
^/Pixi and ^/p2X2 simultaneously, the destination receives 
the symbol y given by 



y 



^Pixi 



P2X2 + Z, where z ^ CSCQ (0,cr^ 



(17) 

and 9i, 62 are i.i.d. random variables distributed uniformly 
over (0, 27r). We assume that only the destination has the 
knowledge of Oi and 6*2. If User-i has the knowledge of 6i, 
then the phase offset can be compensated by each user which 
in turn results in the signal model discussed in Section HI] 
The CC capacity region for the above model can be computed 



A. Numerical Results 

For the equal average power case (i.e., SNR = Pi = P2 
with (T^ = 1), we compute the values of a* for QPSK, 8- 
PSK and 8-QAM constellations. The corresponding values of 
a* are listed in Table |III] The CC sum capacities for QPSK 
and 8-PSK are also provided in Fig. [5] and Fig.|6] respectively 
for the following two cases: (/) the CPA scheme with a = 
a* , and (;7) the CR scheme. For the CR scheme, though the 
two users employ rotated constellations, the channel induced 
rotations (due to 9i and 62) make the effective angle of rotation 
a random variable. Therefore, the CR scheme for this channel 
corresponds to the CPA scheme with a — 1. The figures show 
that the CPA scheme provides larger CC sum capacities than 
the CR scheme at high SNR values. 

For the unequal average power case, we have computed 
a* for the QPSK constellation and for several values of 
SNR = 4 (where ct^ = i). We use the following three 
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NiN2 



fcl=0 fc2=0 



/ {^/P^Xl + y^2X2 : y I 01, ^2) = log2(^iiV2)- 



exp (— |zP/ct^) 



Wi-l N2-1 

E E 

.fci=0 ^2=0 



Afi-lAf2-l 

^ exp(-|/ii(fci,fc2,«i,*2)lV2'T^) 

Lii=0 12=0 



/ii(fci,fc2,ji,*2) = e''i (v/^xi(fci)- /F^a;i(ii)) + e^'^^ (/P^a;2(fc2) - /P^a;2(j2)) 



(18) 
(19) 

(20) 



relations between Pi and P2: (/) P2 = 0.5Pi, (//) P2 = 0.75Pi 
and (//;) P2 = 0.9Pi. The corresponding CC sum capacities are 
presented in Fig. |7]for the CPA and the CR scheme. From the 
figures, note that the CPA scheme provides increased CC sum 
capacities at moderate to high SNR values, especially when 
Pi and P2 are close. Further, the advantage of CPA is noticed 
to diminish when the power difference is more than 3 dB. 



4 




5 10 15 20 25 30 

SNRindB 



Fig. 5. CC sum capacity of QPSK with equal average power constraint and 
random phase offsets. 




5 10 15 20 25 30 

SNRindB 



Fig. 6. CC sum capacity of 8-PSK with equal average power constraint and 
random phase offsets. 



Pg - 0.5P, = 0.75P, 




10 20 30 10 20 30 

SNRindB SNRindB 

Pj = OOP, 




Fig. 7. CC sum capacities of QPSK with unequal average power constraint 
and random phase offsets. 



VI. Conclusion and Directions for Future Work 

The proposed CPA scheme can be useful for two-user 
fading MAC especially when the channel state information 
is available at the transmitters. For such a case, depending 
on the instantaneous channel states of the two users, we 
obtain one of the channel models discussed in this paper 
Therefore, whenever the fading amplitudes of the two channels 
are approximately close, the CPA scheme will be beneficial. 
Hence, studying the impact of the CPA scheme for fading 
MAC is an interesting direction for future work. On the other 
hand, the CPA scheme can also be studied for GMAC with 
arbitrary number of users. However, large number of users in a 
NO-MA scheme is known to increase the decoding complexity 
at the destination [9]. To avoid this increased complexity, the 
proposed results for two-user GMAC can be incorporated by 
the network schedulers, wherein instead of separating all the 
users, the scheduler can separate pairs of users, and then 
employ the CPA scheme between the two users in each pair. 
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